Negative refraction and left-handed behavior in two-dimensional photonic crystals 
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We systematically examine the conditions of obtaining left-handed (LH) behavior in photonic 
crystals. Detailed studies of the phase and group velocities as well as the phase Up and group Ug 
refractive index are given. The existence of negative refraction does not guarantee the existence 
of negative index of refraction and so LH behavior. A wedge type of experiment is suggested that 
can unambiguously distiguinsh between cases of negative refraction that occur when left-handed 
behavior is present, from cases that show negative refraction without LH behavior. 
PACS numbers: 78.20.Ci, 41.20.Jb,42.25.-p, 42.30.-d 



Recently, there have been many studies about, 
rials that have a negative index of refraction rm 
These materials, theoretically discussed by VeselagoH, 
have simultaneoiiS|-|negative permittivity e and perme- 
ability /i. PendryQu has studied structures that they can 
have both e and ^ negative and suggested that they can 
be used to fabricate a perfect lens. In these materials, fc, 
E, H form a left-handed (LH) set of vectors (i.e., S ■ k < 
0, where S is the Poynting vector) and, therefore, are 
called left-handed materials (LHM). Futhermore Markos 
and SoukoulisEI used the transfer matrix technique to cal- 
culate the transmission and reflection properties oLieft- 
handed materials. These scattering data were usedEJ to 
determine the effective permitivity and permeaiulity of 
the LHM and the index of refraction was foundli3 to be 
unambiguously negative in the frequency region where 
both of the retfkwed, values of e and are negative. 
All experimentad'B'0'Qtl that showed left-handed behavior 
were performed in the microwave or millimeter regime. 

In addition Notomi0 studied light propagation in 
strongly modulated two dimensional (2D) photonic crys- 
tals (PCs). Such a photonic crystal behaves as a mate- 
rial having an effective refractive index UeS controllable 
by the band structure. In these PC structures the per- 
mittivity is periodically modulated in space and is posi- 
tive. The permeability /i is equal to one. Negative rics 
for a frequency range was Jound. The existence of neg- 
ative n was demonstratedtll by a finite difference time 
domain (FDTD) simulation. Negative refraction on the 
interface of a 3D PC structure has been experimentally 
observed by Kosaka et al.E3. A negative refractive index 
associated to the negative refraction was reported. Also 
in Ref. 12, large beam steering was observed that was 
called the "superprism phenonoepon" . Similar unusual 
light propagation was observedllj in ID and 2D diffrac- 
tion gratings. In addition, in Luo et. al. it was shown 
that the phenomenon of negativojefraction in the PC can 
be utilized to make a superlensEil. Finally a theoretical 
worlol3 predicted a negative refraction index in PCs. 

Both in the PC and LHM literature there is a lot of 
confusion about what are the correct definitions of the 
phase and group refractive index and what is their rela- 



tion to negative refraction. It is important to examine if 
left-handed behavior can be observed in photonic crys- 
tals at optical frequencies. In addition, it is instructive 
to see how the latter relates with the sign of the phase 
and group refractive index for the PC system. 

In this paper we will attempt to clarify the differences 
between phase refractive index rip and group refractive 
index rig in PC structures. We show that left-handed be- 
havior in photonic crystals is possible and a comparison 
of this LH behavior to regular LH materials is given. We 
will also show that a frequency range exists for which neg- 
ative refraction can occur while it does not correspond to 
regular LH behavior. A wedge type of i^xperiment is sug- 
gested that can unveil the sign of S ■ (i.e. the "right- 
ness" ) in the PC and so the correct sign for the "effective" 
index (phase index Up as we will see in the following). 
This way we can unambiguously distinguish between the 
case of "negative refraction" that corresponds to LH be- 
havior in the PC structure {rip and S ■ Ml3 are negative), 
from cases that show "negative refraction" without LH 
behavior {rip and S ■ k are positive). FDTDEZI simulations 
are presented which support our claims that left-handed 
behavior can be seen in the the PC only in the cases 
where the predicted value for rip (for the infinite system) 
is negative. In all FDTD simulations we used perfectly 
matched layer (PML)ll3 boundary conditions. 

Our results are general but for compjleteness and for 
comparison with the results of Notomitil we studied a 
2D hexagonal lattice of dielectric rods with dielectric con- 
stant 12.96 (e.g., GaAs or Si at 1.55 microns) for the H 
(TE) polarization. The radius of the dielectric rods is 
r — 0.35a, where a is the lattice constant. In Fig. la the 
band structure of the system is shown. In the figure as 
well as in the following the quantity / is the frequency 
normalized as uja/2'Kc, where c is the velocity of light. To 
study the system the surfaces (contours in the 2D system) 
of equal frequency in fc-space,- called equifrequency sur- 
faces (EFS)-, are needed. The EFS basically consists of 
the allowed propagation modes within the crystal for a 
certain frequency. For both the band structur(e-|pad the 
EFS plots, the plane wave expansion methodHc^ with 
1003 plane waves was used. The EFS can have different 
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shapes for various frequencies. We notice that in some 
cases their shape can be almost circular (i.e., uj ~ Ak, 
where A a constant). Basically, these cases are for fre- 
quencies close to the band edges and we will refer to their 
dispersion relation Lu{k) as normal dispersion from now 
on. 

Before we proceed in presenting the FDTD simulation 
results it is useful to discuss the appropriate expressions 
for the phase {Vp) and group (vg) velocities and their as- 
sociated phase (rip) and group (ng) refractive indices for 
the infinite photonic crystal system. The analysis that 
follows is not used in the interpretation of the FDTD re- 
sults, that we discuss later in this paper. However it does 
serve as a theoretical prediction to be compared with the 
FDTD data. For any general case Vp = ^f-j-^ , with 

k = k/k . Also It/gl^lVfeCjI^i^. The value of \np\ for 
a certain angle of incidence 9 will be c\kf{d)\/u! (the re- 
fracted wave- vector kf can be determined from the EFS 
as it will be discussed later) . The sign of Up is determined 
from the behavior of EFS. If the equifrequency contours 
move outwards with increasing frequency then Vg-k >0; if 
they move inwards Vg-k <0. It can be proven analytically 
that for the infinite PC systcpipUjic group velocity coin- 
cides with the energy velocityOu. Therefore the sign of 
Vg ■ k is equivalent to the sign of S ■ Ml3 - where k is in tha 
1"* BZ-. So, to be in accordance with the LH literatureO 
the sign of Up will be the sign oivg-k. In the special case 
of normal dispersion, Up is independent of the angle of 
incidence . For the cases with dispersion close to normal 
Up vs. / for up to the 6*'* band has been calculated and 
the results can be seen in Fig. lb. We observe in Fig. lb 
that \np\ can have values less than one (i.e., superluminal 
phase velocity). However \ng\ has to be always greater 
than one, because the group velocity is the energy veloc- 
ity in the PC, which is less than c. We have calculated the 
group velocity with the k ■ p perturbation methodE3 and 
checked that indeed \ng\ >1 . For the cases with normal 
dispersion Vg—^k with ng=u}d\np\/ duj + From 
this relation and from the graph of Up versus frequency 
(see Fig. lb), it can be checked that indeed, ng >0 when 
Up >0 and Ug <0 when Up <0. This means that for the 
normal dispersion cases negative refraction occurs when 
the phase index is negative. We note this finding is dif- 
ferent of that stated in Ref. 26. For the anisotropic dis- 
persion cases Vg and k are not along the same direction. 
In these cases we choose the sign of Ug to manifests the 
sign of refraction at the interface and therefore Ug will 
have the sign of {vg ■ ki)x {x being the direction along the 
interface and ki the incident wavevector). We would like 
to stress that Ug cannot be used in a Snell-like formula to 
obtain the signal's propagation direction in the crystal. 
However, Up can be used in Snell's law to determine the 
signal's propagation angle, in the normal dispersion case, 
under certain conditionsEll. 

There are two characteristic cases for which an EM 
wave refracts in the "wrong way" (negatively) when it 



hits the PC interface. One can predict the direction of 
the refracted signal from the equifrequency contours in 
the air and PC media. This rne±hology has been used 
widely in ID diffraction gratingscj. The EFS for two fre- 
quencies, -/ — 0.58 and / = 0.535 - that are representa- 
tive for these two cases, are seen in Fig. 2a and Fig. 2b. 
The k/ / conservation condition given by the dashed ver- 
tical line in the figures determines the allowed refracted 
wavevector kf. In both cases there are two choices for 
kf (to point either towards A or towards B). For every 
choosen fcy, the corrresponding Vg will be perpendicular 
to the EFS at that point and point towards increasing 
values of frequency. The correct choice for kf will be the 
one that gives a Vg that points away from the source. 
After checking the EFS for a region around / = 0.58 and 
/ = 0.535 we determine that the frequency increases in- 
wards for the case of Fig. 2a (i.e. Vg ■ k < Q and Up <0 
) while it increases outwards for the case of Fig. 2b (i.e. 
Vg-k > and Hp > 0). Taking this into consideration the 
correct choice for fc/ is to point towards A for the case 
of Fig. 2a and towards B for the case of Fig. 2b. For an 
infinite PC the group and the energy velocity are equal, 
so the group velocity vector (orange vector) in the fig- 
ures represents the direction of propagation for the EM 
signal in the PC. We have confirmed the theoretical pre- 
diction of negative refraction with a FDTD simulation 
for both cases. In other words, we have determined both 
thcorcticaly and numerically that the group index ng for 
both the cases is negativcED. However this type of FDTD 
simulation cannot decide what is the "rightness" of the 
system. 

Therefore we consider an incoming EM wave incident 
normally on a wedged structure, as was the case of the 
UCSD experimentu. We shopki mention that couplings 
to higher order Bragg wave£3 when the wave hits the 
wedged interface are unavoidable in some cases. How- 
ever, the sign of the angle of the 0*'' order outgoing Bragg 
wave coincides with the sign oi-,S ■ k inside the wedge, 
and determines the "rightness"l3 (positive for the right 
and negative for the left handed case) and so the sign of 
Hp of the PC. 

In particular, we show the results for this type of sim- 
ulation for / ~ 0.58 (case of Fig. 2a) and / ~ 0.535 
(case of Fig. 2b) in Figs 3b and 4b respectively. For each 
of these cases the symmetry directions of the fiijst and 
the wedged interface are choosen appropriately^. Nat- 
urally, the k// wavevector will be conserved across both 
interfaces. In Fig. 3b the final outgoing beam is in the 
negative hemisphere (negative in respect to what is ex- 
pected for a regular homogeneous material with positive 
index). This means that the perpendicular to the first 
interface component of the wavevector has reversed sign 
when crossing the air PC interface. Since S is always 
pointing away from the source S ■ k will be negative in 
the PC for the case of Fig. 3b. We note that the second 
outgoing beam in Fig. 3b is due to consequent multi- 



3 



reflection on the upper corner (one can determine that 
by looking at the field inside the wedge). To extract the 
sign of iS • fc in Fig. 4b, because of-the presence of the 
2"'^ beam (higher order Bragg wavecS) we need the addi- 
tional knowledge of the theoretically expected magnitude 
of the angle for the 0*^ order Bragg wave. We use the 
EFS contours for the infinite system to determine the 
magnitude of the refracted wavevector and from that the 
magnitude of the 0*'' order Bragg wave. (We obtained a 
value of 57°). Although the actual magnitude of the re- 
fracted wavevector may differ from the predicted one, no 
dramatic difference is expected. We then observe in the 
FDTD simulation of Fig. 4b that the beam with deflec- 
tion angle of magnitude close to 57° lies in the positive 
hemisphere. So, in the case of Fig. 4b 5 • fc > 0. 

The ray schematics in Figs 3a and 4a - that accompany 
the simulations of Figs. 3b and Fig. 4b- show the theo- 
retically predicted refracted and final outgoing beams by 
using the properties of the infinite system. The beams 
are determined in a similar manner as in the ray schemat- 
ics of Fig. 2. In the following we compare the theoretical 
expected angles (magnitude and sign) with the ones mea- 
sured in the FDTD simulation. The expected outgoing 
angle for the case of Fig. 3b is —36", while the FDTD one 
is ^ —30". Likewise for the case of Fig. 4 the formula in 
Ref. 28 predicts an outgoing wave at 57° and a 1"** order 
Bragg wave at —14°. These values are close to those of 
~ 60° and ^ —10° measured in the FDTD simulation. 

By looking on other bands we can find cases with 
Up > and Ug > 0, as well as rip < and rig > , 
so all combinations of signs for the group and phase in- 
dex are possible. By performing FDTD simulations (of 
both types (slab and wedge)) we have shown that in all 
cases, that the "rightness" of the PC - i.e. the sign of 



coincides with the the theoretical prediction for 
the sign of Tip , clS defined earlier in the paper for the 
infinite PC-, and is irrelevant to the sign of Ug. For the 
case of Fig. 2a, the PC is left-handed and the perpen- 
dicular component of the wavevector reverses sign in the 
air PC interface. In addition it has close to normal dis- 
persion and a negative effective index. So one could say 
that the negative refraction in this case is similar to the 
one occurring ip,fiL.njoub|ixJ(H material that has both e 
and fi negatiyemimWBn 

In conclusion, we have systematically examined the 
conditions of obtaining left-handed behavior in photonic 
crystals. We have demonstrated that the existence of 
negative refraction is neither a prerequisite nor guaran- 
tees a negative effective index and thus LH behavior. Al- 
though we have focused on the presentation of 2D pho- 
tonic crystals, the results are general enough and are ex- 
pected to be correct for 3D PCs. Special care must be 
taken for the couplings of the TE modes and TM modes 
in the 3D systems. We have suggested experiments to 
check our theoretical predictions both in the microwave 
regime as well as in the optical regime. The present work 
will help in the understanding and design of the optical 
devices. 
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FIG. 1: a) The band structure for the H (TE) polarization for 
the PC structure under study, b) The effective phase index 
Up plotted for the frequency regions where the dispersion is 
close to normal. 




FIG. 2: Schematics of refraction at the PC interface for the 
four different cases. The initial (ki) and refracted (fc/) wave 
vectors as well as the group velocity (vg) of the refracted wave 
are drawn. The dashed line represents the conservation of the 
parallel component of the wave vector. 




FIG. 3: Refraction of monochromatic EM wave with / = 0.58 
incident normally on a PC crystal wedge. In (a) the expected 
ray-paths are shown. The wave vectors of the incident (ki) as 
well as the refracted (fc/) and the outgoing (fc//) EM waves 
are shown. In (b) the actual FDTD simulation with the inci- 
dent and outgoing beam is shown. The bold solid line is the 
normal to the wedged interface. 




FIG. 4: The same as figure 3 but for / = 0.535 



